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Abstract 

We study the Lorenz model from the viewpoint of its accessible singulari- 
ties and local index. 
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1 Introduction 

The Lorenz model 



has been studied by Tabor and Weiss (1981) in detail [T5] . In this paper, we study 
the phase space of (CQ) from the viewpoint of its accessible singularities and local 
index. 

As the conditions which each accessible singular point can be resolved, we obtain 
the following 



(1) 



' dx 
~dt 
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y — aex, 



xy — ebz 



—xz + x — ey 



(2) 



e(b - 1)(& - 2a) (b + 3a - 1) = 0, 
(6-1)(6-3<t+1) = 0, 

{b 2 - 5b - 2 - 3(6 - 2)a}{e 2 (b - 1)(76 - 15(7 + 2) - 9} = 0, 
e(b-2a)(b + 3a - 1) = 0. 



These equations can be solved as follows: 



{(a, e, b) = Q, e, o) , (1, -3, 2), (1, 3, 2), (2, 0, 1)}. 
In each case, we study its first integrals, general solutions and phase space. 



2 Accessible singularities 

Let us review the notion of accessible singularity. Let B be a connected open domain 
in C and 7r : W — > B a smooth proper holomorphic map. We assume that H C W 
is a normal crossing divisor which is flat over B. Let us consider a rational vector 
field v on W satisfying the condition 

v G H°(W, e w (-logW)(W)). 

Fixing to £ B and P G Wt , we can take a local coordinate system ( ) of 

W to centered at P such that 7i sm0 oth can be defined by the local equation x± — 0. 
Since -u G iJ°(W, W (- log H)(H)), we can write down the vector field v near 
P = (0, . . . , 0, t ) as follows: 



d 



d 



v = m + ai d Xl 1 Xl dx 2 



+ 



a 2 d 



+ • • • + 



x\ dx r 



This vector field defines the following system of differential equations 



(3) 



dxi 
~dt 

d,X2 



a 1 (x 1 ,x 2 , ....,£„,£), 
a 2 (xi,x 2 , ....,x n ,t) 



Xi 



dx n _ a n (xi,x 2 , ....,x n ,t) 
, dt x\ 

Here a,i(xi, . . . , x n , t), i — 1, 2, . . . , n, are holomorphic functions defined near P = 
(0,...,0,t ). 

Definition 2.1. With the above notation, assume that the rational vector field v on 
W satisfies the condition 



(A) deH°(W,e w (-\ogH)(H)). 
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We say that v has an accessible singularity at P = (0, . . . , 0, to) if 

X\ = and 0,(0, ... , 0, to) = for every i, 2 < % < n. 

If P G 7i smoot h is not an accessible singularity, all solutions of the ordinary 
differential equation passing through P are vertical solutions, that is, the solutions 
are contained in the fiber Wt over t = t . If P G 7Y smoo th is an accessible singularity, 
there may be a solution of (J3]) which passes through P and goes into the interior 

w - n of w. 

Here we review the notion of local index. Let v be an algebraic vector field with 
an accessible singular point p = (0, . . . , 0) and ( coordinate system 

in a neighborhood centered at p . Assume that the system associated with v near 
~p can be written as 



(4) 
d 



dt 



Q 



(x x \ 



%1 



{Q 



a i 



«2 



Q-'-Q 



( Xl\ 



+ 



( xifi(x 1 ,x 2 , 

f2i.X 1 ,X 2} ■ 



■ j ^nj t)^ 



}• 



\ fn\ x li x 2i ■ ■ ■ j x ni t) J 



(f l eC(t)[x 1 ,...,x n ], QeGL(n,C(t)), * G C(t)) 



where fi is a polynomial which vanishes at ~p and fc, i = 2,3, ... ,n are polyno- 
mials of order at least 2 in x\, x 2 , ■ ■ ■ , x n . We call ordered set of the eigenvalues 
) local index at p . 
We remark that we are interested in the case with local index 



(5) 



(l,a 2 /a 1 , . . . ,a n /ai) G Z n . 



If each component of (1, a 2 /ai, a n /ax) has the same sign, we may resolve the 
accessible singularity by blowing-up finitely many times. However, when different 
signs appear, we may need to both blow up and blow down. 

In order to consider the phase spaces for the system ([T]), let us take the com- 
pactification [z : Z\ : z 2 : z 3 ] G P 3 of (x, y, z) G C 3 with the natural embedding 

(,x,y,z) = (z 1 /z ,z 2 /z ,z 3 /z ). 

Moreover, we denote the boundary divisor in P 3 by 7i. Extend the regular vector 
field on C 3 to a rational vector field v on P 3 . It is easy to see that P 3 is covered by 
four copies of C 3 : 

U = C 3 9 (x,y,z), 

I) ) (•V,.^.^) (j = 1,2,3), 
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via the following rational transformations 



X 1 = l/x, Y 1 = y/x, Z 1 = z/x, 
X 2 = x/y, Y 2 = l/y, Z 2 = z/y, 
X 3 = x/z, Y 3 = y/z, Z 3 = l/z. 

The following Lemma shows that this rational vector field v has five accessible 
singular points on the boundary divisor Ti C P 3 . 

Lemma 2.1. The rational vector field v has five accessible singular points: 

'P 1 = {(X 1 ,Y 1 , Z 1 )\X 1 =Y 1 = Z 1 = 0}, 
P 2 = {(X 2 , Y 2 , Z 2 )\X 2 = Y 2 = Z 2 = 0}, 

(6) \ P 3 = {(X 3 ,Y 3 ,Z 3 )\X 3 =Y 3 = Z 3 = 0}, 

P 4 = {(X 2 ,Y 2 ,Z 2 )\X 2 = Y 2 = 0, z 2 = V^T}, 
^P 5 = {(X 2 ,Y 2 ,Z 2 )\X 2 = Y 2 = 0, Z 2 = -^l}. 



Next let us calculate its local index at the points P« 



1,2,3). 



Singular point 


Type of local index 


Pi 


(0,^-1,-^-1) 


Pi 


(0,0,0) 




(0,0,0) 



We see that there are no solutions which pass through P i; (i = 1, 2, 3), respectively. 

In order to do analysis for the accessible singularities P4, P5, we need to replace 
a suitable coordinate system because each point has multiplicity of order 2. 

At first, let us do the Painleve test. To find the leading order behaviour of a 
singularity at t — ti one sets 

X oc 



y oc 



z oc 
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(t 
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(t 


-ti)"' 




c 


(t 


-ti) p ' 



from which it is easily deduced that 

m — 1, 

and 



n 



±2V-1, 



n 



V 



-2. 
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Each order of pole (m, n, p) suggests a suitable coordinate system to do analysis for 
the accessible singularities -P4, -P5, which is explicitly given by 



{X,Y,Z) 



1 y_ z_ 

X x z x z 



In this coordinate, the singular points P4, P5 are given as follows: 



P A = {{X,Y,Z) 

Pr = {(X,Y,Z) 



1-1 1 

~2~' 2 



}• 



— -)> 
2 '2'' 



Next let us calculate its local index at each point. 



Singular point 


Type of local index 


P 4 


(-^1,-2^1,-^1) 


P 5 


(^,2v^T, v^T) 



Now, we try to resolve the accessible singular points -P4, -P5. 
Step 0: We take the coordinate system centered at P4: 



p = X, q = Y 



1 



r = Z — 



r. 



2 1 2 
Step 1: We make the linear transformation: 

p^ 1 ' = p, = q — y/—lr, A 1 ' 

In this coordinate, the system ([1]) is rewritten as follows: 

n \ /jiW 

-2y/=l 


(4,2). This property suggests that we will blow up four times to the direction 
and two times to the direction 

Step 2: We blow up at the point P 4 = {(p {1 \ q (1 \ r (1) ) = (0,0,0)}: 






By considering the ratio ( 1 



+ ...}. 



'1,4,2), we obtain the resonances 



P 



(2) 



P 



(1) „(2) _ £ 



9 



P 



(i) 

or 



.(2) 



r (i) 
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Step 3: We blow up at the point {(p( 2 \ q^ 2 \ r^) = (0, f (6- 1), y/^le(b-2a)}: 

v m- v m J ( 3 )_ g (2) -|(fc-i) r( s) _ ^ ~ V^<b - 2a) 

P P > * p {2) ' p (2) 

Step 4: We blow up along the curve {(p^ 3 \ q^ 3 \ r^)\p^ = 0, 

g(3) = VEI (g2( 6 _ 1)( 76 _ l 5a + 2 ) _ 9)}; 

p (4) =p ( S ) ) w _ g (3) -^(^-l)(7fe-15a + 2)-9) > r(4)=r(3) _ 

' "(3) 



Step 5: We blow up along the curve {{p^ A \ q^\ r^)\p^ = 0, 
= | e (6 _ i) g (4) _ + 2 ) ( £ 2 (6 - i)(76 - 15a + 2) - 9)}: 



u = p^ , V 



- (|e(6 - l)rW - £e(6 + 2) (e 2 (6 - 1)(76 - 15a + 2) - 9)) 



p(4) 



(4) 

w — r ( ' . 



In this coordinate, the system (JTJ) is rewritten as follows: 
(7) 

f du , N 
— = (/i (u, «,«;), 

4 v / ^l£ 3 (6- 1)(6-2(t)(6 + 3a-l) 2e 2 18(6- l)(6-3a + l)w 



dt 9 u 2 27 u 

_ 2e 2 {b 2 - 5& - 2 - 3(6 - 2)a}{e 2 (6 - 1)(76 - 15a + 2) - 9} 

~ ~27 u 
dw y^le 2 (b- 2a)(6 + 3a - 1) 



+ (7 2 (u,u,«;), 



>. rfi 3 u 



+ g 3 (u,v,w) 



where gi(u, v, w) G C[it, t> , w] (i = 1, 2, 3). 

Each right-hand side of the system ((7|) is a polynomial if and only if 



e(b - 1)(6 - 2a)(6 + 3a - 1) = 0, 
(6-1)(6-3<t+1) = 0, 

{6 2 - 56 - 2 - 3(6 - 2)a}{e 2 (6 - 1)(76 - 15a + 2) - 9} = 0, 
e(6-2a)(6 + 3a- 1) = 0. 



These equations can be solved as follows: 

(9) {(a, e, 6) = Q, e, Oj, (1, -3, 2), (1, 3, 2), (2, 0, 1)}. 
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3 The case of (cr, e, 6) = Q, e, 0) 



f dx 
~dt 



(10) 



V x, 

y 3 > 



dt 



— = — xz + x — ey, 
at 



xy. 



This system is equivalent to the 3-rd order ordinary differential equation: 

2 



d 3 x e 3 2 dx Ae f dx\ ldxd 2 x Aed 2 x 
dt 3 3 dt 3x \dt J x dt dt 2 3 dt 2 



This system does not appear in the Chazy polynomial class. 



Theorem 3.1. The phase space X for the system ( JTOl) is obtained by gluing three 
copies of C 3 : 

U j <*C*3{{x j ,y i ,z j )}, j = 0,1, 2 
via tne following birational transformations: 



(12) 

0) x =x, 2/o =3/, = 2, 



1) x 1 =-, 2/i 



(y- 3^ 



lz + K — J -)x + — (3z + e 2 - 3) ) x, 

9 9 



zx — z (3x — 4v^— Te)x, 

6 



2) x 2 =-, 2/2 

x 



(y - + v 73 !^ 



-l(5e 2 + 9), 4e, 2 , , 
+ - (3z + e 2 - 3) ) x, 



z 2 = 2 — (3x + 4v^— Te)x. 
6 



These transition functions satisfy the condition: 

dxi A dyi A dzj = dx A dy A dz (i — 1, 2). 

Theorem 3.2. Let us consider a system of first order ordinary differential equations 
in the polynomial class: 

dx , . dy dz . . 

-^ = h\x,y,z), —=f 2 (x,y,z), — = f 3 {x,y,z). 



We assume that 

(Al) deg(fi) = 2 with respect to x, y, z. 

(A2) The right-hand side of this system becomes again a polynomial in each 
coordinate system (xi,yi,Zi) (i = 1,2). 
Then such a system coincides with the system ( flOl) . 



4 Modified Lorenz model 

In this section, we present 4-parameter family of modified Lorenz model explicitly 
given by 



(13) 



doc s I 

= -gar + y + — (8(9«i - ea 3 ) + V /Z l(24a 2 + a 2 - 16e 2 )) , 

^ = -xz - (24a 2 + - 8£ 2 ) z - + ^(a 3 - A\^le)z 

+ -^(24a 2 a 3 + « 3 - 432«i£ + 64a 3 e 2 - 2 V /Z le(120a 2 + 5a 2 . - 16e 2 ) 

$ = xy + (l2(6ai - ea 3 ) + v /Z T(24a 2 + a 2 )) x-^(a 3 - AV^le)y 



432 



(«3 - 4 v /z Te)(12(6a 1 - a 3 e) + v /Z T(24a 2 + a 2 )), 



where aij, e are complex parameters. 

Theorem 4.1. TTie phase space X for the system ( TT3l) obtained by gluing three 
copies ofC 3 : 

I ) = ^ > {(.<>//;.-;)}• j =0,1, 2 

ma the following birational transformations: 
(14) 

0) x =x, y = y, z = z, 



1) aci =-, |/i 



(w X 

vy 3 



2i = 2 - -(3x - a 3 )x, 
6 



2) x 2 



— 4e 
-lz + ati)x + -gfi z + a 2) ) x, 



1 



x 



y 2 = _((y - -x + V-lz + — (36«i + 24V-la 2 + V-la 3 - 48v-le )x 



36 



4 £ i 
+ —(32 + -(3a 2 + 2 V /Z la 3 £ + 8e 2 ))z, 



z 2 = z — -(3s — a 3 + 8v— le)x. 
6 



These transition functions satisfy the condition: 

dxi A dyi A dzi = dx A dy A dz (i — 1, 2). 

Theorem 4.2. Let us consider a system of first order ordinary differential equations 
in the polynomial class: 

dx dy dz 

= h{x,y,z), — = f 2 {x,y,z), — = f 3 (x,y,z). 

We assume that 

{Al) deg(fi) = 2 with respect to x,y,z. 

(A2) The right-hand side of this system becomes again a polynomial in each 
coordinate system (xi,yi,Zi) (i = 1,2). 
Then such a system coincides with the system ffT3]) . 



5 The case of (a, e, b) = (2, 0, 1 



(15) 



dx 

-dt =m 

dy 

— = —xz + X, 
dt 

dz 



dt 



xy. 



Proposition 5.1. This system has 

(16) I:=x 2 -2z 

as its first integral. 

This system can be solved by reduction to 2nd-order ordinary differential equa- 
tion 

(17) 

or equivalent ly 
(18) 




d 2 x 

~d¥ 



1 3 / 1 I ^ 

~2 +( + 2 )X ' 



--x 6 + (1 + -)x, 



2 v 2' 
which is a special case of equation Ince-Vlll. 
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6 The case of (a, e, b) = (1, 3, 2) 



(19) 



eft 
dy 

dt 

dz 

~dl 



y ~ 3x, 

—xz + x — 3y, 

xy — 6z. 



Proposition 6.1. This system has 

(20) / := e 6t (x 2 - 2z) 

as its first integral. 

Owing to its first integral, this system can be reduced to the system 

dx 

(21) 




(/ + 2e 6t ) 



By making the change of variables 
(22) 



x = y—x, 



Y = 



2x + 4-\x 2 - 2y 

2x ' 



we obtain 



(23) 



^ = X 2 - XY - 2X, 
dt 

dY T 
— = Y 2 - 3XY -2Y- -e- 6t . 
dt 2 



Proposition 6.2. After a series of explicit blowing-ups at ten points including seven 
infinitely near points, the phase space X for this system can be obtained by gluing 
four copies of C 2 x C: 

Uj x C C 2 x C 3 {(xj, yj, t)}, j = 0, 1, 2, 3 
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(*) intersection number 
Figure 1: The bold lines denote (— 2)-curve. 



via the following birational transformations: 
0) x =x, y = y, 



(24) 



1) x l = -, y l 

x 



2) x 2 =-, 2/2 

X 



3) x 3 =xy, y 3 



(yx + ^e~ 6t )x + ^e" 6 * 



((y - 2x)x - -e &t )x + ^-e 



x, 
I 



-at 



X. 



y 



Here, for notational convenience, we have renamed X, Y to x, y. 

We remark that the phase space X is not a rational surface of type E^p (see 
Figure 1). 

It is still an open question whether integrability status of this system is known 
or not. 



7 The case of (a, e, b) = (1, -3, 2) 



(25) 



( dx 

~dl 

dy 

dt 

dz 

~dt 



y + 3x, 

—xz + x + 3y, 
xy + 6z 



Proposition 7.1. This system has 

(26) / := e~ m (x 2 - 2z) 

as its first integral. 
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